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Dipole excitations below the neutron threshold in neutron rich Sn isotopes are studied theoretically 
in the Quasiparticle-Phonon Model with Hartree-Fock-Bogoliubov single particle input. Of special 
interest are the low-lying two-phonon 1~ states and the Pygmy Dipole Resonance (FDR). The 
evolution of low-energy dipole excitations with neutron excess is investigated over the Sn isotopic 
chain including the experimentally unknown regions close to ^'^^Sn. A dependence of the FDR 
strengths and centroid energies on the neutron skin thickness is found. Despite significant multi- 
phonon contributions to mean energies and transition strengths, the FDR states retain their one- 
phonon character. The fragmentation pattern is reduced with increasing neutron excess towards 
the N=82 shell closure which will be of advantage for future experimental work. 

PACS numbers: 21.60.Jz, 24.30.Cz, 27.60.+j 



A genuine feature of neutron-rich nuclei is the appear- 
ance of low-energy electric dipole strength, seen recently 
in high-precision photon scattering experiments already 
in stable nuclei with small 0, and moderate 01 neu- 
tron excess. These so-called Pygmy Dipole Resonances 
(PDR) are observed as a clustering of states close to 
the neutron threshold which in presently accessible stable 
medium- and heavy-mass nuclei is at excitation energies 
~ 5.5—8 MeV. In a first attempt for experimental 
investigations on the dipole strength distribution in ex- 
otic nuclei with radioactive beams recently performed in 
MSU and GSI the neutron-rich (stable and unstable) oxy- 
gen isotopes are studied. A dipole strength is observed 
below Q and above the neutron threshold. The mea- 
surements of the latest in the energy region 5—8 MeV 
for i8'20q )g^j,g consistent with that for other nuclei that 
may exhibit the Pygmy Dipole Resonance'. Although 
carrying only a small fraction of the full dipole strength 
the PDR states are of particular interest because they 
are expected to reflect the motion of the neutron skin 
against the core of normal nuclear matter. In order to 
obtain more direct evidence for such a peculiar mode 
studies of low-energy dipole excitations over sufficiently 
long isotopic chains are necessary. Here, we present re- 
sults of an exploratory theoretical investigation for the 
neutron-rich Sn isotopes in experimentally less or even 
unknown regions at neutron numbers N=70 — 82. As dis- 
cussed in (3i , the nature of the low-energy dipole strength 
differs significantly from the isovector Giant Dipole Reso- 
nance (GDR) mode where proton and neutron fluids as a 
whole move against each other. Moreover, the PDR mode 
has to be distinguished from the other known low-energy 
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isoscalar dipole excitation, namely the two-phonon 1^ 
states resulting from the anharmonic interactions of the 
lowest 2"*" and 3~ states in a nucleus. The anharmonici- 
ties are reflecting the intrinsic fermionic structure of the 
nuclear phonons thus deviating from ideal bosons. These 
competing effects are taken into account in an appropri- 
ate way by the Quasiparticle Phonon Model (QPM)[^. 
Applications of the QPM to low-energy dipole strength 
0iA .9i IQJ E^nd also in the recent PDR investigations in 
^°®Pb Q have led to very good descriptions of data thus 
giving confidence on the reliability of the model for such 
investigations. 

The QPM approach is at present the only method al- 
lowing for sufficiently large configuration space such that 
a unified description of low-energy single and multiple 
phonon states and the GDR is feasible. In the one- 
phonon sector the Quasiparticle Random-Phase Approx- 
imation (QRPA) is recovered as a well understood limit- 
ing case if the quasi-boson approximation is used. Such 
a unified treatment is exactly what is required in order 
to separate the multi-phonon and the genuine PDR 1~ 
strengths in a meaningful way. For the aim of this paper 
we use the standard form of QPM i.e. approximat- 
ing residual interactions in terms of separable multipole- 
multipole interactions with empirical coupling constants, 
see e.g. iHi, 12]. The isoscalar quadrupole and octupole 
coupling constants are chosen to reproduce the exper- 
imental energies and electromagnetic transitions of the 
2^ and 3j~ states. Since for A<126 where the data are 
available the coupling constants vary smoothly with neu- 
tron number on a level of less than 5% extrapolations into 
the unexplored mass region towards ^^^Sn can be done 
safely. 

Since single particle energies and a reliable description 
of ground state properties in general are critical quanti- 
ties for extrapolations of QRPA and QPM calculations 
into unknown mass regions here we put special emphasis 
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on the mean-field part. Because of the numerical con- 
straints set by the QPM a semi-microscopic approach 
is chosen. Taking advantage of the Kohn-Sham theo- 
rem [T^ of Density Functional Theory (DFT) the total 
binding energy B(A) is expressed as an integral over an 
energy density functional with (quantal) kinetic (t) and 
self-energy parts, respectively, 



B{A) ^ J 



(1) 



and 
The 



pairing contributions are indicated by j^pair- 
second relation is obtained from Hartree-Fock- 
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Bogoliubov (HFB) theory with occupancies Vj and po- 
tential energies < U >j of the occupied levels j, see e.g. 
[T^ITsl . Above, U{p) is the proper self-energy, i.e. not in- 
cluding the rearrangement contributions from the intrin- 
sic density dependence of nuclear interactions |0, flGj . 
Hence, U{p) has to be distinguished from the effective 
self-energy obtained by variation 



1 dpUjp) 

2 dp 



(2) 



and appearing in the single particle Schroedinger equa- 
tion. In order to keep the QPM calculations feasible we 
choose S = Uws to be of Wood-Saxon (WS) shape with 
adjustable parameters. By inversion and observing that 
the densities and potentials in a finite nucleus are nat- 
urally given parametrically as functions of the radius r, 
we find 



p{r)U{r) ^ -2 



'ds^Uwsis) . (3) 
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Evaluating these relations with the microscopic pro- 
ton and neutron densities obtained by solving the 
Schroedinger equation with Uws the potential U{p) is 
the self-consistently derived reduced self-energy entering 
e.g. into the binding energy. 

In practice, for a given nucleus of mass A the depth 
of the central and spin-orbit potentials, radius and dif- 
fusivity parameters of Uws are adjusted separately for 
protons and neutrons to the corresponding single 
tide separation energies, the total binding energy 
the charge radii and (relative) differences of proton and 
neutron root-mean-square (RMS) radii, 



6r = \/ 



(4) 



from our previous HFB calculations [11 m. The theo- 
retically obtained RMS radii are compared to those deter- 
mined from charge exchange reactions by Krasznahorkay 
et al. [isl for a number of Sn isotopes. However, 
since the measurements did not provide absolute 6r val- 
ues the data were normalized in [isL IT^ to theoretical 
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FIG. 1: Ground state properties of the Sn isotopes. In the 
upper panel the total nuclear binding energies per particle cal- 
culated with the adjusted WS potentials are compared to the 
values from the Audi-Wapstra compilation p^. In the lower 
panel the calculated with the adjusted WS potentials differ- 
ences of proton and neutron rms radii are compared to the 
experimental values obtained by Krasznahorkay et al. from 
charge exchange reactions in ref . iljjj . 



predictions. Hence, strictly spoken only the relative mass 
dependence of the experimental results shown in Figgis 
of significance. The HFB energies were scaled by an av- 
erage effective mass of m* /m = 0.68 thus removing the 
known problem of unrealistically large HFB level spac- 
ings at the Fermi surface. The approach sketched above 
leads to very satisfactory results on binding energies and 
proton-neutron RMS-differences as shown in Fig^for the 
Sn isotopes. A smooth dependence of the parameters on 
A is found which supports the reliability of the method. 

In the QPM calculations the structure of the excited 
states is described by a wave function including up to 
three-phonon configurations Q which are built from a 
basis of QRPA states of natural parity excitations with 
J'^ = 1~ 5~. Since the one-phonon configurations up 
to E:r=20 MeV are considered the core polarization con- 
tributions to the transitions of the low- lying 1~ states 
are taken into account explicitly. Hence, we do not 
need to introduce effective charges. The two- and three- 
phonon configurations are truncated at Ea;=4.5 MeV for 
the calculation of the quadrupole-octupole 1~ state for 
the proper comparison with the available Nuclear Reso- 
nance Fluorescence (NRF) data |7j]. For the QPM calcu- 
lations between 4.5-^8 MeV the two- and three-phonon 
basis is limited to states up to E2.=8.5 MeV and 8 MeV, 
respectively. 

As an important prerequisite for the reliability of the 
model wave functions we consider first the lowest J'^=2+ 
and 3'^=3- states in the i20^i30gj^ isotopes. The QPM 
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TABLE I: QPM results for the energies and the reduced 
B(E1), B(E2) and B(E3) transition probabihties of the first 
1", 2"*" and 3" states in ^^""^'"'Sn isotopes. A comparison 
with the experimental data 7] is presented. 
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energies and transition probabilities of these states are 
presented in Table I and compared to the known exper- 
imental data. The well reproduced properties of the 2^ 
and 3j~ states are important for a trustable description of 
the two-phonon quadrupole-octupole J'^^lj" state which 
is of our particular interest. The experimental data on 
excitation energies and transition strengths l2(||, avail- 
able for ^^^^^•^''Sn, are well described by the QPM as 
seen from Table L 

The QPM predicts for the lowest-lying 1^ states 
excitation energies varying from Ea;=3.203 MeV to 
E^=4.115 MeV in izo^iaogj^ ,^ggg n^^^^e I). In these iso- 
topes the 1^ states are at a level of about 90% predom- 
inantly given by two-phonon quadrupole-octupole con- 
figurations. Since with increasing neutron number the 
energy of the latter becomes larger, certain other, higher- 
lying, two-phonon configurations also start to contribute. 
The three-phonon configurations are most important for 
the lower-mass Sn isotopes because of their open shell 
structure. The reduction correlates with the decrease 
of collectivity when approaching the N=82 shell closure. 
The QPM results for the boson forbidden El transitions 
l^j from the two-phonon 1~ state to the ground state and 



between the 3^ and 2^ states are also shown in Table |T] 
The experimental data available for ^^^Sn-^^^^'Sn are 
well described. 

Of central interest for this work are the 1~ states 
above the two-phonon dipole state and below the neu- 
tron threshold and their evolution with the neutron ex- 
cess. These states are eventually to be identified as PDR 
modes related to excitations of the neutron skin as found 
in ^"^Pb 3]. In fact, in all nuclei considered here the first 
three 1" QRPA states contribute to this energy region. 
They are rather well separated by a energy gap of more 
than 1.3 MeV from other higher- lying one-phonon 1~ 
states which are more likely to belong to the low-energy 
tail of the GDR. QRPA and QPM results on the total 
PDR strengths and average centroid energies (defined by 
the ratios of energy and non-energy weighted sum rules) 
in Sn isotopes are compared in TabUll As a general re- 
sult we find a correlation of the total PDR strength and 
the neutron excess as refiected by the increase of the to- 
tal PDR strength with mass number. The PDR centroid 
energies, however, decrease with mass number. They are 
only weakly affected by multi-phonon admixtures. 

The calculated energies and transition st reng ths are in 
a good agreement with experimental data fl(||. The de- 
tailed analysis of the El strength distribution over the 
1~ excited states in -^^°^^'^°Sn isotopes reveals that the 
major part of the strength is concentrated in the states 
with relatively large one-phonon [1^] or [1;^] components. 
For the total PDR strengths the differences between the 
QRPA and QPM are significant (see Table|lll). Although 
the contribution of the higher-lying one-phonon compo- 
nents to the wave function structure of the low- lying 1^ 
states is small their infiuence on the total El transition 
is important because of partial admixtures of the large 
collective one-phonon matrix elements from the GDR re- 
gion. The QPM B(E1) values are enhanced significantly 
by factors ranging from about 2 to 1.5 for A — 120 
to A — 130, respectively. Thus, the PDR transition 
strengths contain a considerable fraction of strength from 
higher-lying (GDR) states, virtually admixed through 
the anharmonic QPM interactions. These admixtures, in 
fact, overcompensate the loss of ground state transition 
strength resulting from the dissipation of the one-phonon 
states into the multi-phonon components. In ^^^Sn and 
^^°Sn the QPM calculations exhaust about 85% of the 
PDR one-phonon strength in the bound energy region. 
In Table HTl also the relative contribution of the PDR to 
the Energy Weighted Sum Rule (EWSR) is given. 

An interesting question is to relate the PDR modes to 
the neutron skin. While in 3J this was done on the level 
of transition densities here we analyze the PDR strength 
distributions as functions of <5r, eq. 2| We find two well 
separated regions with a sudden increase in the slope be- 
yond A = 126 where 5r starts to exceed 0.2 fm. Hence, 
a clear signal on the thickness of the neutron skin is ob- 
tained allowing to extract information on nuclear shapes 
from PDR data. 

The QRPA and QPM dipole strength distributions be- 
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TABLE II: Reduced transition probabilities and mean energies of the PDR modes in i^OSn-i^'^Sn 

calculated within QRPA 

are presented. The contribution of the PDR mode in the Energy Weighted Sum Rule (EWSR) is given in %. For Sn the 
calculations are compared with the experimental data 



Nucleus (E) [MeV] J2 BjEl) j [eVm^] EWSR (%) 

Exp. QRPA QPM Exp. QRPA QPM 

^^°Sn ~- 6^839 6^920 - 0142 0289" .185 

i^^Sn - 6.695 6.708 - 0.181 0.347 .228 

i2*Sn 6.5'^'" 6.544 6.641 0.345(43)'^°' 0.219 0.398 .267 

i^^Sn - 6.330 6.442 - 0.261 0.416 .304 

"**Sn - 6.100 6.210 - 0.281 0.434 .313 

"°Sn - 5.940 6.150 - 0.337 0.495 .361 

"^Sn - 5.900 - - 0.360 - .380 



low neutron threshold in ^^^Sn and ^^°Sn are compared 
in FigEl The lowest QPM 1" states do not have QRPA 
counterparts because they are the quadrupole-octupole 
two-phonon states. The QPM anharmonicities introduce 
a considerable fragmentation - in these cases up to about 
80 states. Compared to ^3 we find less fragmentation 
because our two- and three-phonon configuration spaces 
are somewhat smaller. However, the obtained for PDR 
total strength and centroid energy is in very reasonable 
agreement with the available experimental data in i24sn 
(see Table inil. 

The lowest QRPA 1~ states are mainly given by pure 
neutron two-quasiparticle excitations with small admix- 
tures of higher-lying two-quasiparticle neutron and pro- 
ton configurations on a level of a few percent or less. The 
dominance of neutron excitations in this energy region is 
in agreement with the results of previous PDR studies 
in Sn and other nuclei by relativistic random-phase ap- 
proximation (RRPA) nil, DFT [13 and QPM 0. The 
remarkable stability of the wave function structure over 
the isotopic chain resembles what is found theoretically 
also for the GDR. However, we emphasize that the PDR 
states are of a much stronger isoscalar (or eventually 
mixed isospin) content than the GDR, see also |^]. Hence, 
the PDR states cannot be considered as simply being the 
low-energy tail of the isovector GDR, despite the anhar- 
monic admixtures discussed above. Rather, these dipole 
modes are of a genuine character which cannot be de- 
duced by extrapolations from the GDR region. This 
conclusion is supported by the differences in transition 
densities and velocity fields discussed in Q . 

In summary, neutron rich tin isotopes were studied in 
a semi-phenomenological approach combining HFB and 
QPM theory. The results confirm the success of such 
a microscopically inspired description. An important 
step in understanding the dipole spectra is to disentangle 
the PDR states from the low-energy two-phonon dipole 
states, achieved here by using the QPM approach with 



up to three-phonon configurations. For ^^°Sn ^^^Sn 
we obtained low-energy dipole strength concentrated in 
a narrow energy interval such that a pygmy dipole res- 
onance (PDR) can be identified. The evolution of the 
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FIG. 2: QRPA and QPM calculations of the El strength dis- 
tribution in a/i^^Sn and b/^^^Sn below the neutron threshold. 



PDR strength distribution with neutron excess shows 
that transition strengths and energy locations are indeed 
closely correlated with the neutron skin. Hence, mea- 
surements of PDR strength distributions will provide in- 
formation on nuclear shapes far off stability. 
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